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Supermassive black holes (SMBHs) and their host galaxies are generally thought 
to coevolve, so that the SMBH achieves up to about 0.2 to 0.5% of the host 
galaxy mass in the present day. The radiation emitted from the growing SMBH 
is expected to affect star formation throughout the host galaxy. The relevance 
of this scenario at early cosmic epochs is not yet established. We present spec¬ 
troscopic observations of a galaxy at redshift z = 3.328, which hosts an actively 


1 


accreting, extremely massive BH, in its final stages of growth. The SMBH mass 
is roughly one-tenth the mass of the entire host galaxy, suggesting that it has 
grown much more efficiently than the host, contrary to models of synchro¬ 
nized coevolution. The host galaxy is forming stars at an intense rate, despite 
the presence of a SMBH-driven gas outflow. 

Several lines of observational evidence, spanning a wide range of cosmic epochs, have led to 
a commonly accepted picture wherein supermassive black holes (SMBHs, Mbh > 10® Mq; Mq 
is the solar mass) coevolve with their host galaxies (i^). Moreover, energy- and/or momentum- 
driven “feedback” from accreting SMBHs (Active Galactic Nuclei; AGN) is thought to quench 
star formation in the host galaxy (5). To directly test the relevance of such scenarios at early 
cosmic epochs (high redshifts, z) requires the most basic properties of SMBHs and their hosts, 
including masses and growth rates, to be observed. Several observational studies found that 
at z ^ 2 (more than 3.3 billion years after the Big Bang), the typical BH-to-stellar mass ratio, 
Mbh/A/*, increases towards higher redshifts (6-8), suggesting that some SMBHs were able 
to gather mass more efficiently, or faster, than the stellar populations in their hosts. To date, 
measurements of Mbh at earlier epochs (z > 2) have only been conducted for small sam¬ 
ples of extremely luminous objects [Tagn > 10"^®ergs“^ (9-12)] representing a rare subset 
of all accreting SMBHs, with number densities of order 1 to 10 per Gpc® [i.e., ~ 10 ® to 
10“® Mpc“® (/3)]. Moreover, the high AGN luminosities in such sources overwhelm the host 
galaxy emission and prohibit a reliable determination of M*, and therefore of Mbh/M*. We ini¬ 
tiated an observational campaign aimed at estimating Mbh in x-ray-selected, unobscured z ^ 3 
to 4 AGN within the Cosmic Evolution Survey field [COSMOS; (14)]. Such sources have lower 
AGN luminosities and are more abundant than the aforementioned luminous sources by factors 
of 100 to 1000 (13,15) and thus form a more representative subset of the general AGN pop¬ 
ulation. Moreover, the fainter AGN luminosities and rich multiwavelength coverage of AGN 
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within the COSMOS field enable reliable measurements of the mass and growth rate of the 
stellar populations in the host galaxies (M* and star-formation rate, SFR). 

CID-947 is an x-ray-selected, unobscured AGN at z = 3.328, detected in hoi]\XMM-Newton 
and Chandra x-ray imaging data of the COSMOS field [see Fig. [S4] and sections IS^ and IS4l in 
the supplementary materials {16)]. We obtained a near-infrared (IR) ^-band spectrum of CID- 
947 using the MOSFIRE instrument at the W. M. Keck telescope, which dX z = 3.328 covers 
the hydrogen H/9 broad emission line (see details in section[ST|in the supplementary materials). 
The calibrated spectrum shows a very broad H/9 emission line, among other features (Fig.[T]). 
Our spectral analysis indicates that the monochromatic AGN luminosity at rest-frame 5100 A 
is T 5100 = 3.58lo;o8 x 10"^® ergs“k The typical line-of-sight velocity, i.e. the full-width at half¬ 
maximum of the line, is llddOlgyo kms“^ (see section [SI . 21 in the supplementary materials). 
By combining this line width with the observed T 5100 and relying on an empirically calibrated 
estimator for Mbh, based on the virial motion of ionized gas near the SMBH (77), we obtain 
^BH = 6 . 9 I 5 2 X 10® Mq. All the reported measurement-related uncertainties are derived by 
a series of simulations and represent the 16th and 84th quantiles of the resulting distributions. 
These simulations indicate a SMBH mass larger than 3.6 x 10® Mq at the 99% confidence level 
(see sections IS 1 .2! and IS^ for more details). Determinations of Mbh from single-epoch spectra 
of the H/7 emission line are known to also be affected by significant systematic uncertainties, of 
up to ~ 0.3 to 0.4 dex. For a detailed discussion of some of the systematics and related issues, 
see [JS3]in the supplementary materials. This high Mbh is comparable with some of the most 
massive BHs known to date in the local universe (75),or with the masses of the biggest BHs in 
the much rarer, more luminous AGN at 2 ; ~ 2 to 4 [e.g., (9)]. The bolometric luminosity of CID- 
947 is in the range Lboi — (1.1 — 2.2) x 10"^® erg estimated either from the observed optical 
luminosity or the multiwavelength spectral energy distribution. Combined with the measured 
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Mbh, we derive a normalized accretion rate of L/LEdd — 0.01 to 0.02. This value is lower, by at 
least an order of magnitude, than the accretion rates of known SMBHs at z ^ 3.5 [e.g., (9, 10)]. 
Further assuming a standard radiative efficiency of 10%, we obtain an e-folding time scale for 
the SMBH mass of at least 2.1 x 10® (Gy; see section IS^ . which is longer than the age of the 
universe at z = 3.328. By contrast, even the most extreme models for the emergence of “seed” 
BHs predict masses no larger than M^eed ~ 10® Mq at z ~ 10 to 20 [e.g., (19)]. Therefore, the 
SMBH powering CID-947 had to grow at much higher accretion rates and at a high duty cycle 
in the past, to account for the high observed Mbh only 1.7 Gyr after z ~ 20. CID-947 could 
have evolved from a parent population similar to the fast-growing SMBHs observed in 2 ; ^ 5 
quasars, which have L/LEdd ~ 0.5 to 1 and Mbh — 10® Mq [e.g., (11,12)]. The requirement 
for a high accretion rate in the very recent past is supported by the clear presence of a high- 
velocity outflow of ionized gas, observed in the rest-frame ultraviolet spectrum of the source 
(fig.[S4l). The broad absorption features of C IV A1549 and Si IV A1400 have maximal velocities 
of i^max — 12,000kms“\ Assuming that this outflow is driven by radiation pressure, these 
velocities require accretion rates of L/LEdd ~ 0.1, as recently as 10® to 10® years before the 
observed epoch (see section [S4l) . We conclude that the SMBH powering CID-947 is in the final 
stages of growth and that we are witnessing the shut-down of accretion onto one of the most 
massive BHs known to date. 

The rich collection of ancillary COSMOS multiwavelength data available for CID-947 en¬ 
ables us to study the basic properties of its host galaxy (see details in section [S^ in the supple¬ 
mentary materials). A previously published analysis of the observed spectral energy distribution 
of the emission from the source reveals an appreciable stellar emission component, originating 
from 5.6loi x 10^® Mq in stars (20). Our own analysis provides a yet lower stellar mass, of 
M* = X 10^® Mq. However, we focus on the previously determined, higher stellar 

mass, as a conservative estimate. The source is also detected at far-IR and (sub)millimeter 
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wavelengths, whieh allows us to eonstrain the SFR in the host galaxy to about 400 Mq year“^ 
The stellar mass of the host galaxy is eonsistent with the typieal value for star-forming galaxies 
at 2 ; ~ 3 to 4 [i.e., the “break” in the mass funetion of galaxies; (27)]. Similarly, the eombina- 
tion of M* and SFR is eonsistent with the typieal values observed at 2 ; ~ 3 to 4, whieh appear 
to follow the so-ealled main sequenee of star-forming galaxies (22). Thus, the host galaxy 
of CID-947 is a typieal star-forming galaxy for its redshift, representing a population with a 
number density of about 5 x 10“® Mpc“^ [e.g., (27)]. This suggests that neither the intense, 
ionizing radiation that emerged during the fast SMBH growth, nor the AGN-driven outflow, 
have quenehed star formation in the host galaxy. The relatively high stellar mass and SFR of 
the host galaxy further suggest that it is unlikely that the AGN affeeted the host in yet earlier 
epoehs. That is, even in this ease of extreme SMBH growth, there is no sign of AGN-driven 
suppression of star formation in the host. 

Our analysis indieates that the BH-to-stellar mass ratio for CID-947 is Mbh/A 7* — 1/8. In 
eomparison, most loeal (dormant) high-mass BHs typieally have Mbh/A 7* ~ 1/700 to 1/500 
[see Fig. [ 2 ] and, e.g., (4, 23)]. The Mbr/M^ value that we find for CID-947 is thus far higher 
than typieally observed in high-mass systems in the loeal universe, by at least an order of mag¬ 
nitude and more probably by a faetor of about 50. The only loeal system with a eomparably 
extreme mass ratio is the galaxy NGC 1277, whieh was reported to have Mrr/M^ ~ 1/7 [with 
Mbh = 1-7 X 10^° Mq ~ 2.5 X Mbh(CID-947); see (24), but also (25)]. At earlier epoehs 
(still z < 2), the general trend is for Mbh/M* to inerease slightly with redshift, but typieally 
not beyond Mbh/M* ~ 1/100 (see Fig.[3|). Only a few systems with reliable estimates of Mbh 
show Mbh/M* reaehing as high as 1/30 [e.g., (<5-S)]. 

Given the high masses of both the SMBH and stellar population in CID-947, we expeet 
this system to retain an extreme Mbh/M* throughout its evolution, from z = 3.328 to the 
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present-day universe. Beeause the Mbh that we find is already eomparable to the most massive 
BHs known, it is unlikely that the SMBH will experienee any further appreciable growth (i.e., 
beyond Mbh — 10^° Mq). Indeed, if the SMBH accretes at the observed rate through z = 2, 
it will reach the extreme value of ~ 10^° Mq, and hy z = 1 it will have a final mass of ~ 
2.5 X 10^° Mq. As for the host galaxy, we can constrain its subsequent growth following several 
different assumptions. First, if one simply assumes that the galaxy will become as massive as 
the most massive galaxies in the local universe [M* ~ 10^^ Mq', (2(5)], then the implied final 
mass ratio is on the order of Mbh/M* ~ 1/100. Alternatively, we consider more realistic 
scenarios for the future growth of the stellar population, relying on the observed mass (M*) and 
growth rate (SFR). Our calculations involve different scenarios for the decay of star formation 
in the galaxy (see section [S5] in the supplementary materials), and predict final stellar masses 
in the range M* {z = 0) ~ (2 — 7) x 10^^ Mq, which is about an order of magnitude higher 
than the observed mass at z = 3.328. The inferred final mass ratio is Mbh/M* ~ 1/50. This 
growth can only occur if star formation continues for a relatively long period ( ^ 1 Gy) and at 
a high rate (> 50 Mq year“^). This would require the presence of a substantial reservoir, or the 
accretion, of cold gas, which, however, could not increase the SMBH mass by much. Finally, 
in the most extreme scenario, the star formation shuts down almost immediately (i.e., due to 
the AGN-driven outflow), and the system remains “frozen” at Mbh/M* ~ 1/10 throughout 
cosmic time. If the SMBH does indeed grow further (i.e., beyond 10^° Mq), this would imply 
yet higher Mbh/M*. Thus, the inferred final BH-to-stellar mass ratio for CID-947 is, in the 
most extreme scenarios, about Mbh/M* ~ 1/100, and probably much higher (see Fig.©. 

CID-947 therefore represents a progenitor of the most extreme, high-mass systems in the 
local universe, like NGC 1277. Such systems are not detected in large numbers, perhaps due to 
observational selection biases. The above considerations indicate that the local relics of systems 
like CID-947 are galaxies with at least M* ~ 5 x 10^^ Mq. Such systems are predominantly 
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quiescent (i.e., with low star-formation rates, SFR 1 Mq year and relatively rare in the 
local universe, with typical number densities on the order of 10 ^ Mpc ^ (2(5). We con¬ 
clude that CID-947 provides direct evidence that at least some of the most massive BHs, with 
Mbh ~ 10^° Mq, already in place just 2 Gy after the Big Bang, did not shut down star forma¬ 
tion in their host galaxies. The host galaxies may experience appreciable mass growth in later 
epochs, without much further black hole growth, resulting in very high stellar masses but still 
relatively high Mbh/M*. Lower-mass systems may follow markedly different coevolutionary 
paths. However, systems with Mbh/M* as high as in CID-947 may be not as rare as previously 
thought, as they can be consistently observed among populations with number densities on the 
order of ~ 10“^ Mpc“^, both at z > 3 and in the local universe, and not just among the rarest, 
most luminous quasars. 
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Figure 1 : The observed Keck/MOSFIRE spectrum and best-fit model for the H/? emission 
complex of CID-947. The data are modeled with a linear eontinuum (dotted), a broadened 
iron template (dot-dashed) and a eombination of broad and narrow Gaussians (dashed), whieh 
eorrespond to the H/) and [O ill] emission lines (see seetion lST^ in the supplementary materials 
for details regarding the speetral modeling). The broad component of H/3 has a full width at 
half maximum of FWHM (H/)) = 11330 kms“^, which results in Mbh = 6.9 x 10® Mq and 
Mbh/^* = 1/8. The red dashed line illustrates an alternative scenario, in which the SMBH 
mass derived from the H/) line width would result in Mbh/M^, = 1/100 [i.e., FWHM (H/)) = 
3218 kms“^], clearly at odds with the data. The spike at Arest — 4640 A is due to a sky feature. 
The bottom panel shows the residuals of the best fit model. 
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galaxy stellar mass, (Mq) 

Figure 2: A comparison of CID-947 with a compilation of observed Mbh and M* esti¬ 
mates in the local universe [adapted from (4), assuming the tabulated bulge-to-total fraetions]. 
CID-947 (red star) has a very high BH-to-stellar mass ratio of Mbh/M^ ~ 1/10. The asym- 
metrie error bars shown on Mbh and M* represent measurement-related uneertainties, while 
the symmetrie ones demonstrate systematie uneertainties of 0.3 dex (on Mbh) and 0.1 dex (on 
M*). The masses inferred for subsequent growth seenarios are highlighted as empty red stars. 
The CID-947 system is expeeted to evolve only mildly in Mbh (perhaps to 10^° Mo), but 
M* should grow to at least 2 x 10^^ Mq, and possibly to as mueh as ~ 7 x 10^^ Mq, by 2 ; = 0. 
The loeal galaxies NGC 1277 and M87, whieh eould be eonsidered as deseendants of systems 
like CID-947, are highlighted as filled symbols [ (25) and (27), respeetively]. Some studies 
suggest these galaxies to have somewhat higher Mbh, and therefore relatively high mass ratios, 
of Mbh/M* = 1/7 and 1/127, respectively (24,28). 
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Figure 3: The observed cosmic evolution of the BH-to-stellar mass ratio, Mbh/ 3^*? and 
its extrapolation beyond z ^ 2. CID-947 (red star) has M^h/M^ = 1/8 at 2 ; ~ 3.3, 
whieh is higher by a faetor of at least ~ 50 than the typieal value in loeal, inaetive galax¬ 
ies (at most, Mbh/ 3F* ~ 1/500; dotted line). The error bars shown for CID-947 represent 
only the measurement-related uneertainties, propagating the uneertainties on Mbh and on M*. 
The different data points at z < 2 represent typieal (median) values for several samples with 
Mbh/ 3T* estimates, with uncertainties representing the scatter within each sample [filled sym¬ 
bols, open circles, and open triangles represent samples from (7), (29), and (6), respectively; 
adapted from (7)]. Even compared to the extrapolation of the evolutionary trend supported by 
these lower-redshift data, Mbh/M^ ^ {z + 1)^ [dashed line, scaled as in (30)], CID-947 has a 
significantly higher Mbh/3T*. 
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mass (M*) and star formation rate (SFR; ^ jS2l) : the derivation of black hole mass (Mbh) and 
accretion rate (L/LEdd), and a discussion of the SMBH evolution (^|S3]); the rest-frame UV 
spectrum, broad absorption features and properties of the AGN-driven outflow and of the 
calculations of possible final M* of the host galaxy ( ^}S5]) . 

SI New Keck/MOSFIRE K-hand Data and Analysis 

Sl.l Keck/MOSFIRE Observations and Data Reduction 

The source CID-947 (J2000.0 coordinates a = 10:01:11.35, 6 = -1-02:08:55.6) was observed 
with the Keck/MOSFIRE instrument (31) during the night of January 23-24th., 2014, with 
observing time allocated through the Yale-Caltech collaborative agreement. We used the normal 
K-band setup, which covers order 4 of the 110.5 mm“^ reflection grating. As CID-947 was our 
primary target, it was positioned near the center of the mask, providing a spectral coverage 
of A =19,415-23,837 A. To ensure adequate coverage of the sky background emission, and 
its subtraction from the AGN signal, we used 4 pairs of MOSFIRE bars, to form a 24"-long 
pseudo-slit. The MOSFIRE pixel scale is 0.18"/pix. To prevent significant slit losses, we set 
the slit width(s) to 1", which resulted in a spectral resolution of about R = A/AA = 3600. 
Observational conditions during the night were generally good, with typical seeing of ~0.8" in 
the A-band during the science observations. The science exposures, totaling an hour, consisted 
of 20 separate sub-exposures of 3 minutes each, dithered between two positions along the slit 
with a separation of 4", to allow for an accurate subtraction of the sky emission. The typical 
airmass during the observations was about 1.08. Several times during the night we also observed 
the AOv stars HIP-34111 and HIP-56736, as well as the fainter white dwarf GD71, to allow a 
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robust flux calibration. 


The data were redueed using a eombination of different tools. First, we used the dedieated 
MOSFIRE pipeline (2014.06.10 version) to obtain flat-fielded, wavelength ealibrated 2D spee- 
tra of all the sourees observed within eaeh mask (ineluding the standard stars). The wavelength 
ealibration was performed using sky emission lines, and the best-fit solutions aehieved an rms 
of ~ O.lA. Next, we used standard I RAF proeedures to produee a ID speetrum, using an aper¬ 
ture of 11 pix (i.e., 2"). Finally, we used the Spextool IDL paekage to remove the tellurie 
absorption features near 2 /rm and to perform the relative and absolute flux ealibrations, based 
on a detailed library speetrum of Vega (32, 33). The absolute flux ealibration we obtained is in 
exeellent agreement with the arehival photometry available for CID-947: the synthetie magni¬ 
tude derived from the speetrum is = 20.03 (AB magnitudes), eompared with the arehival 
value of 20.00 ± 0.01 (34). We however ehose to apply the minor sealing needed to mateh the 
arehival photometry (a faetor of 1.03), in order to be fully eonsistent with the value used in the 
SED deeomposition (^ }S2l) . We finally note that CID-947 is one of several COSMOS targets 
observed in this eampaign, and we verified the robustness of the different reduetion steps by 
visually verifying that the same reduetion yields artifaet-free speetra for the other sourees. The 
typieal signal-to-noise aeross the eore part of the speetrum is S/N ~ 5 — 7 per pixel. Af¬ 
ter re-binning the speetrum to a uniform spaeing of 1 A (in rest-frame; 60kms“^), we obtain 
S/N ^ 7 — 10 per speetral resolution element. 


S1.2 Spectral Analysis of H/3 Emission Complex 


We modeled the ^-band speetrum of CID-947 to measure the monoehromatie eontinuum lumi¬ 
nosity at (rest-frame) wavelength of 5 lOOA (\Lx 5100A , or L 5100 ) and the width of the broad 
H/3 emission line (FWHM (H/3)). The analysis methodology is very similar to that diseussed in 
numerous previous works (e.g., (35-37), and referenees therein), and is only briefly deseribed 
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here. 


The speetra were modeled with a linear (pseudo) eontinuum, a broadened Fe ii template 
(35), and a eombination of Gaussians whieh aeeount for the broad and narrow emission lines, 
namely H/5, [Oiii] AA4959,5007, and Hell A4686. The H/? model eonsists of a broad eom- 
ponent (modeled with 2 Gaussians) and a narrow eomponent, whieh is tied to the width of 
the [Out] lines. The eontinuum flux at SIOOA was estimated direetly from the best-fit lin¬ 
ear eontinuumQ We preferred to use FWHM over cTblr as the probe of the virial veloeity 
field of the BLR gas, as the former ean be more robustly estimated in speetra of limited 
S/N, as is the ease with our MOSFIRE data (39). However, sinee the best-fitting model for 
H/S turned out to be overwhelmingly dominated by a single broad Gaussian eomponent (the 
flux ratio between the two eomponents is 45:1), the differenees between the results derived 
by two approaehes are expeeted to be negligible. The best-fit models are presented in Fig¬ 
ure |ST] The relevant best-fit parameters resulting from our fitting of the MOSFIRE speetrum 
are FWHM(H/S) = 11330kms“\ and L 5100 = 4.16 x lO'^^ergs”^. After aeeounting for 
host-galaxy eontamination of about 14% (following the analysis presented in ^ jS2l) . the intrinsie 
optieal luminosity beeomes L 5100 = 3.58 x 10"^^ergs“^. We note that the broad H/3 profile 
in our best-fit model is highly symmetrie, with an asymmetry index of A.I. = 0.03, eonsistent 
with the typieal value found in large samples of un-obseured AGN (see, e.g., (40)). The broad 
eomponent is, however, blue-shifted by about 830kms“^, relative to the expeeted wavelength 
(at the souree systemie redshift). Sueh large blue-shifts are relatively rare, with an oeeurrenee 
rate of only about 5%y To verify that our estimate of FWHM (H/5) is not severely affeeted by 
these properties of broad H/3 profile, we obtained an alternative estimate of the line width, whieh 
relies only on the blue part of the profile, whieh is not blended with Iron and [O ill] emission. 

^The monochromatic luminosity, as all other luminosities and ages reported here, are calculated assuming a 
cosmological model with Ha = 0.7, Hm = 0.3, and Hq = 70kms“^ Mpc“^. 

^This estimate is based on the H/3 measurements of about 20,000 SDSS AGN, presented in Ref. (37). 
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This alternative estimate is obtained by doubling the one-sided line width, which was measured 
following the approach described in Ref. {41). This results in a FWHM (H/3)-equivalent of 
9236 kms“^, only 0.09 dex lower than our best-fit value. However, we stress that this alterna¬ 
tive estimate is far less robust than our fiducial, best-fit value, as it relies on the identification of 
the (spectral) pixel with maximal flux density, which is very sensitive to small-scale flux density 
fluctuations in spectra of limited S/N, as in the present case. We conclude that our estimate of 
the width of the broad component of H/9 is not significantly affected by the shape or shift. 

Due to the complexity of the spectral fitting procedure, the best approach to derive the 
measurement-related uncertainties on L 5100 and FWHM (H/5) (and therefore on Mbh; see ^JS3] 
below) is via re-sampling of the data. To this end, we constructed a set of 500 artificial spectra, 
by adding normally-distributed random noise to the observed spectra of CID-947, scaled to 
provide either S/N = 10 (comparable to the noise level in the observed data). An additional 
set of 500 simulations assumed a more conservative noise level oi S/N = 5, to verify that 
our results are not driven by an under-estimation of the noise level in the data. We then re¬ 
fitted each of these artificial spectra, using the same fitting procedure as described above. The 
resulting best-fitting models are illustrated in the top panel of Figure [SH After measuring L 5100 
and FWHM (H/3) for each of these best-fitting models, we obtained an artificial sample of 
FWHM (H/9) values, the cumulative distribution function of which is illustrated in the lower- 
left panel of Fig.[SlJ Fig. [ST] clearly demonstrates that, even under the conservative assumption 
of S/N = 5, about 95% of our simulations resulted in FWHM (H/9) ^ 8600kms“^. This 
conservative lower limit on FWHM (H/9) is lower than the best-fit value by about 0.14 dex. 

The spectral region adjacent to the [O lll] lines may suggest that their profiles may include 
an additional broad component (i.e., a “blue wing”). We therefore performed yet another set of 
simulations, with an alternative version of the fitting procedure that allows for an additional 
broad component for [Olll]. The broad components for the two [Olll] lines were forced 
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to share a common width (in the range 500 — 1400 km s“^) and relative shift (in the range 
—350 —1-150, kms“^). These limits are motivated by the distributions of line widths and shifts 
found for large samples of un-obscured AGN (e.g., (42-44)). The results of the simulations 
indicate that this adjustment to the [O III] profiles does not systematically affect our estimate of 
FWHM (H/9). The median value obtained in the simulations is FWHM (H/?) = 10150 kms“\ 
and 95% of the simulations resulted in FWHM (H/?) > 6740 kms“^. 


S2 Broad-band Spectral Energy Distribution and Estimates 
of Lbob and SFR 

We used the available multiwavelength data for CID-947 to determine the bolometric luminos¬ 
ity (Lboi) of the AGN in CID-947, and the stellar mass (M*) and star formation rate (SFR) of 
the host galaxy. The broad-band spectral energy distribution (SED) for CID-947 includes data 
from a large variety of surveys of the COSMOS field, including data in the X-ray (Chandra and 
XMM-Newton), optical-to-near-IR (Subaru and CFHT), mid-to-far-IR (Spitzer and Herschel)), 
and millimeter (JMCT) regimes. 

The broad-band SED of CID-947 was analyzed in previous studies of COSMOS AGN, 
which derived and reported estimates of Lboi- One such analysis, based on Chandra X-ray 
data (45), yields Tboi.sEO = 1-31 x 10^® ergs“^ (correcting for our adopted cosmology). □ An¬ 
other study, based on XMM-Newton X-ray data (ref. (4(5); XMM-ID 60131), gives Tboi,SED = 
1.81 X 10^®ergs“^ (also cosmology-corrected). The difference between these two values, of 
a factor of 1.6, is mostly due to the markedly different X-ray fluxes reported for CID-947 in 
the Chandra and XMM-Newton surveys of the COSMOS field, which may be due to intrinsic 
source variability. We note that the stellar component in the SED (see below) has a negligible 


^Here we use the luminosity integrated between 40 keV and 1 fjm, to avoid double-counting the re-processed 
mid-IR emission (which would add about 17% to Lboi for CID-947). 
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contribution to these estimates of Lboi (i-e., < 1%). 

To derive the host galaxy properties, the SED of CID-947 was modeled separately for the 
rest-frame UV-optieal-NIR, and for the Mid-IR-to-millimeter regime. 

The rest-frame UV-optieal-NIR part of the SED ineludes the emission from the aeereting 
SMBH, a part of whieh is re-proeessed by a dusty toroidal strueture (“torus”) and re-emitted in 
IR wavelengths, and from stellar population of the host galaxy. The data in this regime eonsists 
of flux measurements in 13 speetral bands (obtained with Subaru, CEHT and Spitzer), ranging 
from ~ 3700 A (CEHT/m*) to 24 pm (Spitzer/MlPS). We first rely on the data aeeumulated 
and analyzed in a previous COSMOS study by Bongiorno et al. (20). Here we mention briefly 
only some of the features of this modeling, and refer the reader to (20) for a detailed diseussion. 
In that study, the data were modeled as the sum of two distinet eomponents, representing the 
emission originating from the AGN and from stars in the host galaxy. The AGN eomponent 
is deseribed by the multiwavelength AGN SED of Riehards et al. (2006; (47)). The stellar 
eomponent was deseribed by a grid of models, produeed by a well-established stellar popular 
synthesis proeedure (48). Eaeh of the templates represents a stellar population with a different 
age (ranging from 50 Myr - 1.88 Gyr) and exponential deeay rate (tsfh = 0.1 — 30 Gyr), and 
further assume a Chabrier initial mass funetion (49). The templates were then subjeeted to both 
nuelear and galaxy-wide dust extinetion (with Eb-y values of up to 1 and 0.5, respeetively). 
The best-fit model for CID-947 in the Bongiorno et al. study (Eig. [S2l left) provides a stellar 
mass of M* = S.STtgis ^ with a reddening of Eb-y = 0.5. The fraetion of the total 

monoehromatie luminosity at 5100 A whieh is eontributed by the stellar eomponent is about 
/host (5100) = 0.14. This host eontamination is taken into aeeount when we estimate the mass 
of the SMBH (see gS3]). 

We have repeated the UV-to-IR SED fitting, using the most up-to-date imaging data avail¬ 
able for CID-947 (UltraVISTA DR2, (34); see Table |ST]), and a slightly modified AGN model. 
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in which the emission from the dusty torus (dominating the mid-IR regime) is separated from 
the intrinsie AGN radiation (dominating the UV-optieal regime). Our analysis resulted in a very 
similar stellar mass to the aforementioned one. The stellar eomponent is represented, as before, 
by the Riehards et al. AGN SED, whieh is however extrapolated as a power-law at Arest > 4000 
A. The IR emission from the dusty torus is represented by eomposites from a dedieated study 
of the IR SEDs of AGN (50). The grid of stellar population models have remained the same as 
in the aforementioned “referenee” SED fit, with the age of the stellar population eapped at the 
age of the Universe at z = 3.328. Our ehoiee of the Bruzual & Chariot models is motivated 
by the faet that they were also used in most studies of star forming galaxies at 2 ; > 3, whieh 
we use here as referenee (e.g., (21,22)). Eor a detailed diseussion of the effeets of alternative 
stellar population models, e.g., (2J) and (22). We have restrieted the eomponents so that the 
UV-optieal part of the SED would be dominated by the AGN, i.e. /agn > /host- This is mo¬ 
tivated by the overall AGN luminosity of CID-947 (see above) and the faet that the rest-frame 
UV speetrum does not show signifieant host eontamination in the deep absorption features (see 
^JSH). We have also explieitly omitted the data below Arest = 1216 A (i.e., the u and B bands), 
as these are expeeted to be affeeted by Eyct absorption by the intergalaetie medium along the 
line of sight to CID-947. The resulting additional best-fit model (Eig. [S^ right) relies on a 
stellar population with an age of 1 Gyr, a stellar mass of M* = 4.37lo/g ^ 10^° Tf©, with a red¬ 
dening of Eb-v = 0.05. In this new fit, the stellar eomponent eontributes /host (5100) = 0.36. 
Next, we re-fitted the data with a restrieted model in whieh the stellar population is kept at 
the oldest reasonable age (1.8 Gyr). Sinee older stellar populations have higher mass-to-light 
ratios, sueh a fit would in prineiple provide a eonservative upper limit on the stellar mass. This 
fit resulted in M* = b.dSloJg x 10^° Mq, and /host (5100) = 0.35. Einally, we have repeated 
this oldest-population fit, this time without the restrietion of /agn > /host- This resulted in 
M* = 5.95 ± 0.21 X 10^° Mq, and /host (5100) = 0.55. We stress however that these latter 
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age-restricted models do not provide the best fits of the data, as the resulting is higher than 
that found for the non-restricted case. We conclude that the best-fit stellar mass that we obtain 
for the host of CID-947 is M* = 4.37 xlO^° Mq. We however choose to base the rest 

of the analysis on the slightly higher mass found in the Bongiorno et al. study, as it represents 
a more conservative choice given the extremely high BH-to-stellar mass ratio we find for CID- 
947, and since it is based on the same SED decomposition code that was used in some of the 
reference studies to which we compare our results (<5). 

The Mid-IR-to-millimeter part of the SED is dominated by emission from the (cold) dusty 
gas in the host galaxy, heated by the star formation activity. Here we rely on Herschel and 
JMCT/AzTEC detections at 500 /xm and 1.1 millimeter, and upper limits at 100, 160, 250 and 
350 /xm, from Herschel (51,52). The data were fit with a grid of dust-emission templates of star¬ 
forming galaxies, covering a representative range of SED shapes (i.e., effective temperatures; 
(55)). The best-fit template for the far-IR and millimeter data implies a star formation rate of 
SFR = 392M0yr“^. We present the Mid-IR-to-millimeter data and models in Eig.[S3l We 
stress that the AGN contribution to the emission in this regime is negligible, as demonstrated 
by the dotted black line in Eig. [S3] (54). The SER estimate relies on low-resolution Herschel 
and JMCT/AzTEC measurements, and therefore may be contaminated (confused) by emission 
originating from neighboring (unrelated) sources. 

S3 Determination of black hole mass and accretion rate, and 
the past evolution of the SMBH 

Using the best-fit values for L 5100 and FWHM (H/3) (see gS1.2l abovel. and applying commonly- 
used virial estimators of Mbh (10,37), we obtain Mbh = 6.91 x 10® Mq. Alternative calibra¬ 
tions of such virial mass estimators do not alter our findings significantly. For example, the cal¬ 
ibration obtained by a recent reverberation mapping study (55) implies Mbh = 5.68 x 10® Mq; 
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that is, smaller than our fiducial measurement by less than 0.1 dex. We use the simulations de¬ 


scribed above ( ^S1.2I) to estimate the measurement uncertainties on Mbh- The best-fitting L 5100 
and FWHM (H/9) for each simulated spectrum were combined to provide a set of artificial es¬ 
timates of Mbh- The cumulative distribution function of these Mbh estimates is presented in 
the lower-right panel of Fig.[STJ For the S/N = 10 simulations we find that the 16% and 84% 
quantiles are at 5.73 x 10® and 7.66 x 10® Mq, respectively, resulting in Icr-equivalent uncer¬ 
tainties on Mbh of —1.18 x 10® and -1-0.75 x 10® Mq. We further find that 95% (99%) of these 
simulations resulted in Mbh > 5.09 x 10® Mq (3.61 x 10® Mq). The corresponding values for 
the S/N = 5 simulations are Mbh > 3.82 x 10® and 2.96 x 10® Mq, respectively. Virial (or 
“single-epoch”), H/9-based estimates of Mbh are also known to be prone to systematic uncer¬ 
tainties, of up to ~ 0.4 dex, due to the reliance on the empirical T^blr — T 5100 relation, and the 
overall normalization of the mass estimators. We stress however that our analysis of CID-947 
should, in principle, suffer less from systematics, compared to similar studies of more luminous 
sources at 2 ; ~ 3.5, as its luminosity of L 5100 = 3.4 x 10"^^ ergs“^ lies within the range covered 
directly by reverberation mapping experiments (55-57). The leading systematic uncertainty in 
the present case is therefore associated with the assumption of a typical “geometrical factor” 
(commonly referred to as /blr), which is of order 0.1 dex (58,59). 

Another source of concern is the possibility that the BLR is observed at a high inclination 
angle, as suggested by the presence of the broad absorption features (i.e., BAL features) in 
the rest-frame UV part of the spectrum of CID-947 (see ^JS4l below). One may suspect that in 
such a case, the measured FWHM (H/3) would systematically over-estimate the typical velocity 
dispersion in the BLR, and thus lead to an overestimated Mbh- We have investigated this issue 
by comparing the distributions of FWHM (Mg ll) in large samples of non-BAL QSOs, and 
those with a BAL feature in the C IV A1549 line (as in CID-947), among sources at 1.6 < 2: < 
1.9, drawn from a large catalog based on the Sloan Digital Sky Survey (36). Although in the 
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current work we used H/9, and not the Mg It A2798 line as our virial estimator, the widths of 
these two lines have been shown to be elosely eorrelated, and they are thought to originate from 
a similar region within the BLR (e.g., (36, 37)). The distributions of FWHM (Mg It) for BAL 
and non-BAL QSOs (1723 and 15370 objeets, respeetively) are very similar in shape, with 
the median FWHM value for BAL QSOs being only slightly higher, by merely 230kms“\ 
Moreover, there is no exeess of BAL QSOs with line widths eomparable to what we estimate 
for CID-947 (i.e., > 10, OOOkms"^). 

Several phenomenologieal studies have raised the possibility that, for some luminous sourees 
with partieularly broad H/S lines (FWHM (H/9) ^ 4000 km s“^), the line profiles might inelude 
a signifieant eontribution from a non-virialized “very broad eomponent”, whieh should not be 
taken into aeeount when estimating Mbh (e.g., (60)). Simply adopting the empirieally derived 
(and perhaps luminosity-dependent) eorreetions suggested in sueh studies (60,61), our estimate 
for Mbh should be sealed down by about 0.2-0.25 dex. In the eontext of our main finding, of an 
extremely high BH-to-stellar mass ratio for CID-947, this would mean Mbh/M* ~ 1/15—1/12 
- still a very high value (see Fig. [3]). We however stress that the eommonly advoeated approaeh 
to singling out sueh peeuliar objeets is based on the unambiguous identifieation of two eom- 
ponents in the broad H/3 emission line (the “eore” and the “very broad eomponent”), as well 
as some line asymmetry. Sinee the H/3 profile in CID-947 does not show sueh a eomplieated 
strueture, the aforementioned empirieal eorreetions should not be applied. 

Finally, the alternative, “one sided” estimate of the line width (see ^SL2I above) would 
translate to a deerease of about 0.18 dex in any “virial” estimate of Mbh. As explained in 
^SL2[ and demonstrated in our simulations, this is not a robust estimate of the line width. We 
note, however, that sueh a deerease in Mbh would have a similar effect on our main result as 
the one discussed above, namely providing Mbh/M* ~ 1 /12. 

We eonelude that the SMBH powering the AGN in CID-947 has a mass of Mbh > 3.6 x 
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10® Mq, at the 99% confidence level, and our best estimate is Mbh = 6.9 x 10® Mq. 

The rest-frame optical spectrum was used to derive yet another estimate of Lboi, by apply¬ 
ing a bolometric correction (i.e., /boi (siOOA^ = Lboi/Tsioo)- This approach is consistent with 
many previous studies of un-obscured AGN, at all redshifts. We used the luminosity-dependent 
prescription described in (37), which in turn relies on the 5-band bolometric corrections pre¬ 
sented in (62), translated to 5100 A assuming a UV-optical SED with oc v (63). For 
CID-947, this results in Lboi, opt = 2.11 x 10"^® erg s“^. This value is highly consistent with the 
XMM-Newton-hd&QA estimate of Tboi, sed (rcf- (4<5); within 0.06 dex), but significantly higher 
than the C/ianJra-based value (ref. (45); by a factor of almost 2). 

The derived values of Lboi and Mbh were combined to provide estimates of the normalized 
accretion rate, in terms of the Eddington luminosity, L/LEdd = Lhoi/ (1-5 x 10 ^®Mbh/Mq) 
(this definition of the “Eddington ratio” assumes a Solar gas-phase metalicity). We obtain 
L/L'&dd = 0.021 for the Lsioo-based estimate of Lboi, or 0.019 and 0.011 for the XMM-Newton 
and Chandm-hased estimates of Lboi, sed, respectively. Given the fact that the estimates of Lboi 
were obtained using very different approaches, and the systematic uncertainties associated with 
the estimation of Lboi (o-g-, the scatter in /boi (siOOA^), we consider these estimates of L/Lsdd 
to be in excellent qualitative agreement: the SMBH in CID-947 is accreting at a rate of at most 
L/Lsdd — 0.02. 

By combining the estimated L/Lsad and a standard radiative efficiency of ?7 = 0.1 (62), 
we obtain an e-folding timescale for the growth of the SMBH, following the expression tbh = 
4 X 10® E/LEdd resulting timescales are about 2.1 and 4 Gyr, for the higher (Lsioo-based) 

and lower (Chandm-hased) estimates of L/L^dd, respectively. In any case, these timescales are 
longer than the age of the Universe at z = 3.328, of about 1.88 Gyr, and than the elapsed time 
since the earliest seed black holes likely formed, 1.7 Gyr (z ~ 20). This very long timescale 
thus suggests that the SMBH in CID-947 had to experience an earlier epoch of much faster 
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growth (i.e., higher accretion rate), and/or that it had to originate from the most massive type of 
seeds. In particular, extrapolating the growth history “backwards” to z = 10, assuming constant 
accretion rate and radiative efficiency of L/L^dd = 0.02 and = 0.1 (respectively), results in 
^BH = 10) ~ 3.7 X 10® Mq. Even with the lowest efficiency within the framework of a 
geometrically-thin, optically thick accretion disk, rj = 0.038 (maximally retrograde spinning 
BH; e.g., (64)), the implied mass is still about 10® Mq. The most extreme BH seed production 
mechanisms rely on different “direct collapse” scenarios, but generally provide very few seeds 
as massive as M^eed ~ 10® Mq (see, e.g., (65) and reviews in (19, 66,67)). Some very recent 
studies speculate that some BH seed masses may be yet higher, perhaps up to 10® Mq, but 
not before z = 10 (68). Considering the accretion rate onto the SMBH, some recent models 
highlight the possibility of yet more efficient accretion, as the disk becomes “slim” and sur¬ 
passes the simplified (spherical) Eddington limit, perhaps reaching M/M^dd ~ 3 (69). We 
note however that such extreme models for BH seed production and accretion may not neces¬ 
sarily be required to explain objects like CID-947. The implied Mbh of CID-947 at 2 ; ~ 5, of 
about 3 X 10® Mq, is consistent with that observed in the population of high-luminosity quasars 
at that epoch (11, 12). Such sources, however, have much higher accretion rates, typically 
L/L^dd ~ 0.5, and can emerge from standard accretion (at relatively high duty cycles and low 
radiative efficiency), from a broad range of BH seed masses, including those of stellar remnants 
(i.e., Mseed ~ 10® Mq). We stress that in any case, CID-947 had to have higher-than-observed 
accretion rate sometime in the past, to account for its high mass. As we show in ^JS4] below, 
we have evidence that this epoch of high accretion rate took place relatively shortly before the 
observed epoch. 
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S4 Rest-frame UV spectrum and BAL features 


An optical spectrum of CID-947 was obtained as part of the zCOSMOS survey (ref. (70); 
zCOSMOS-ID 823936), and we present it in Fig. [Sll h clearly shows very broad and deep 
absorption troughs blue-ward of the Si IV A1400 and C IV A1549 lines, identifying CID-947 as 
a Broad Absorption Line Quasar, a sub-population that comprises about 10-20% of luminous, 
un-obseured AGN (BALQSOs; see, e.g., (71-73), and references therein). Moreover, the ab¬ 
sorption feature blue-ward of the A1 iii A185 7 line suggests that CID-947 may belong to the yet 
rarer sub-class of low-ionization BAL QSOs (i.e., it is a “LoBAL”). For both C IV and Si IV we 
estimate a maximum outflow velocity of |vmax| — 12, 000 kms“^, with the C IV trough proba¬ 
bly slightly broader. While this value of Umax is not uneommon among BALQSOs (73), it is an 
outlier in the Umax — L/LEdd plane: virtually all known BALQSOs with eomparable Umax have 
mueh higher accretion rates, typically L/L-^dd > 0.1 (74). 

A simple model for the launching of such high-velocity outflows (75, 76) yields a value for 
the maximum outflow velocity Umax — 9300 ^ — O.l) ^ kms“\ 

where 7?abs,o.i is the distance of the absorber from the eontinuum source, sealed to 0.1 pc; Mg 
is the SMBH mass, sealed to 10^ Mq, /q i is the fraetion of continuum photons absorbed (or 
seattered) by the outflowing gas, scaled to 10%; and A ^22 is the absorber column density, sealed 
to 10^^ cm“^. Assuming the observed value for Mbh, and also A^h = 10^^ cm“^, i?abs,o.i = 1 
and /o.i = 1, this expression implies that a wind with Umax — 12, 000 kms“^ should have been 
launehed by a SMBH aeereting at L/LEdd ~ 0.1, and probably at rates as high as L/Lsad — 0.7 
(see (76) for diseussion of viable ranges on all parameters). This is significantly higher, by at 
least an order of magnitude, than the observed value of L/LEdd- We also note that, since within 
the framework of this model L/LEdd oc i?abs,o.i, the implied accretion rate can easily reaeh 
L/LEdd 1 if the outflow has reaehed Ipc. 


30 



The high-velocity outflow was launched at a time At ~ -Rabs/^'max prior to the observed 
epoch. Even for a conservative assumption of i?abs ~ Ikpc (e.g., ref. (77) and references 
therein), the implied age of the outflow is about 10^ years. Any alternative, more realistic 
assumption regarding i?abs would imply an even shorter timescale. 

Thus, the relatively high terminal velocity of the observed outflow lends further support to 
the scenario in which CID-947 was accreting at much higher rates, probably as recently as 
~ 10® years before the observed epoch. 

S5 Subsequent evolution of M* 

We estimated the “final” stellar mass of CID-947 (M* {z = 0)) in several ways, all of which 
rely on the observed stellar mass (5.6 x 10^° M©) and star formation rate (~ 4OOM0yr“^). 
Our calculations assume that we are observing the host galaxy of CID-947 near its peak of star 
forming activity, and that the SFR can only decline with time. 

First, we assume an exponential decline in SFR, with typical e-folding timescales in the 
range of r = 1 — 2 Gyr. These short timescales are supported by several observational stud¬ 
ies which constrain the ages of the stellar populations in massive, low-redshift galaxies (78), 
and basically implies that the final mass of the galaxy can be approximated by M* (x = 0) ~ 
M* (z = 3.328) -f SFR x r. These “integrated” masses should be scaled down, by factors of 
about 1.6, to account for the fact that some of the mass is returned back to the galaxy’s gas 
(e.g., through stellar winds, Supemovae explosions etc.). Using the observed M* and SFR for 
CID-947 we obtain final masses of M* ( 2 ; = 0) ~ (2.8 — 5.3) x 10^^ Mq. 

An alternative calculation relies on the scenario in which star-forming galaxies evolve on the 
“main sequence” at all epochs, until they quench, and that the mass functions of such objects, at 
all epochs, can be linked via the so-called “continuity approach” (79). In particular, we assume 
that for star forming galaxies the specific star formation rate, sSFR, evolves as sSFR l/t. 
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and that the probability of quenching depends predominantly on M*. The different parameters 
in the calculations were derived by applying the continuity approach to the observed (evolving) 
stellar mass functions (see ref. (79) for details). These calculations predict that a galaxy with 
properties like those of CID-947 will evolve to reach 7.2 x 10^^ Mq at z = 0, if it never 
quenches. Alternatively, it can reach 2.2 x 10^^ Mq if it quenches at Xquench = 2.51, when 
the (mass-dependent) probability of quenching reaches 75%. Allowing for some additional 
mass growth through galaxy-galaxy mergers would increase the latter estimate of final mass 
by a factor of about 2.6, to 5.75 x 10^^ Mq. The increase due to mergers for the former mass 
estimate is minute, since the probability of experiencing a similar-mass merger decreases with 
increasing mass. In any case, these calculations show that the stellar population in CID-947 
should grow by about an order of magnitude between x ~ 3.3 and x ~ 0. 

We note that the high SFR we measure in CID-947 implies a molecular gas mass of at least 
M (H 2 ) ~ 2 X 10^° Mq, and perhaps as high as 2 x 10^^ Mq {80, 81). This suggests that the 
host galaxy can experience significant growth by consuming this reservoir of cold gas, without 
any additional gas accretion from its surroundings (i.e., the IGM), nor from mergers. 
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Figure SI: Results of the re-sampling simulations for the H/3 emission eomplex of CID- 
947. Top - (A): Results for 500 simulations, with S/N = 10 and “nominal” [Olll] profiles 
(i.e., single-Gaussian). The diagram illustrates the original speetrum (blue), the best-fit linear 
eontinuum (dashed blaek), and the eolleetion of 500 models that fit the simulated speetra (gray 
shaded region). Bottom: Cumulative distribution funetions (CDFs) of the obtained line widths, 
FWHM (H/S) (B) and SMBH masses, Mbh (C). In both panels, the solid lines illustrate the 
results for the S/N = 10 simulation, while the dashed lines eorrespond to a more eonservative 
simulation, with S/N = 5. The (dotted) horizontal lines mark the 5 and 50% (i.e., median) 
levels. We find that 95% of the S/N = 10 simulations resulted in FWHM (H/3) > 9615 kms“^ 
and Mbh > 5.09 x 10® Mq; and that 99% resulted in FWHM (H/9) > 8175 kms“^ and Mbh > 
3.61 X 10® Mq. 
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Figure S2: The UV-to-IR SED of CID-947, based on the available ground-based and Spitzer 
imaging data in the COSMOS field, whieh is used to infer the stellar mass of the host galaxy. 
Left - SED fitting from Bongiorno et al. (2012; (20)). The observed SED (blaek points) is fit by a 
model (black line) that consists of three distinct components: an un-obscured AGN, dominating 
the UV regime (solid blue) and a stellar population (magenta), which contributes a significant 
fraction of the emission around rest-frame wavelength of ~ 1 pm. Right - our own best-fitting 
model, in which the dusty toroidal structure, dominating the mid-IR regime (dashed blue), is 
treated separately from the intrinsic AGN emission. 
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Figure S3: The mid-to-far IR SED of CID-947, based on the available Spitzer data and the low 
resolution Herschel and millimeter-wave observations, whieh is used to infer the star formation 
rate of the host galaxy. The gray lines represent a subset of the far-IR templates of star-forming 
galaxies we used (53), with the best-fit template and the eorresponding SFR highlighted in 
red. The AGN eontamination at sub-millimeter-to-millimeter wavelengths (Arest > 200 /im) 
is negligible, as illustrated by the pure-AGN spectral energy distribution (black dotted line; 
adapted from ref. (54)). 
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Figure S4: Optical spectrum of CID-947, probing the rest-frame UV regime, obtained as part 
of the zCOSMOS survey (70). Dashed vertical lines mark the expected center wavelength of 
the Si IV A1400, C IV A1549, and A1 III A1857 lines. Blue-shifted, broad absorption features are 
clearly detected next to each of these three lines, identifying CID-947 as a low-ionization broad 
absorption line QSO (or “LoBAL QSO”). The broad absorption troughs blue-ward of Si IV and 
C IV reach Umax — 12000 kms“^. The red part of the spectrum (Arest < 1900 A), including the 
C III] A1909 line, is known to be affected by (instrumental) fringing. 
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Table S1: UV-to-IR Spectral Energy Distribution 

Telescope/ 

Instrument 

band 

-^obs 

rriAE 

i ± A 

-^rest 

[/xm] 

log {XLx] 
[erg s” 

) ± A 

Subaru/ 

V 

5449 A 

20.900 

± 

0.007 

0.1259 

45.937 

± 

0.003 

SuprimeCam 

r 

6232 A 

20.803 

± 

0.006 

0.1440 

45.918 

± 

0.002 


/+ 

7621 A 

20.394 

± 

0.004 

0.1761 

45.994 

± 

0.002 



8801 A 

20.588 

± 

0.002 

0.2033 

45.854 

± 

0.001 

VISTA/ 

Y 

1.020 /im 

20.639 

± 

0.003 

0.2357 

45.770 

± 

0.001 

VIRCAM 

J 

1.250 /im 

20.333 

± 

0.003 

0.2888 

45.803 

± 

0.001 


H 

1.650 /im 

20.290 

± 

0.004 

0.3813 

45.700 

± 

0.001 


Ks 

2.154 /im 

20.029 

± 

0.005 

0.4976 

45.689 

± 

0.002 

Spitzerl 

chi 

3.526 /im 

19.734 

± 

0.011 

0.8147 

45.593 

± 

0.004 

IRAC 

ch2 

4.461 /im 

19.589 

± 

0.009 

1.0307 

45.549 

± 

0.003 


ch3 

5.677 /im 

19.349 

± 

0.073 

1.3117 

45.540 

± 

0.027 


ch4 

7.704 /im 

18.765 

± 

0.075 

1.7800 

45.641 

± 

0.028 

MIPS 

24/i 

23.68 /im 

16.868 

± 

0.026 

5.4704 

45.912 

± 

0.010 
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